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Tests of the Kobayashi-Maskawa model [l| for CP vi- 
olation are ongoing. In particular, B- factories are di- 
recting focus towards measurements of the lesser known 
angles of the Cabibbo-Kobayashi-Maskawa (CKM) tri- 
angle, fa and fa. Measurements of fa typically rely 
on time-dependent CP-violation studies of B meson de- 
cays to 7r + 7r~, p ± ir T and p + p~ since the leading 
tree amplitudes for these processes involve the relevant 
CKM phases. However, penguin amplitudes may also 
contribute significantly in these decays and — via intro- 
ducing additional unknown phases — greatly impair fa 
constraints from the time-dependent measurements. In 
such cases, isospin analyses can be employed to separate 
the tree- level process from penguin contamination 

Measurements of fa from the pn system rely on knowl- 
edge of the B° — ► p°tt° branching fraction 0,0]. Since 
the tree amplitude of B° — > p°7r° decay is color sup- 
pressed, the decay rate is sensitive to the penguin am- 
plitude contribution. Thus, the branching fraction 
plays a critical role in constraining the fa uncertainty due 
to penguin pollution from time-dependent B° — ► p^w^ 
measurements 0, Furthermore, measurement of fa 
from the full B — * pn isospin analysis requires the p°tt° 
branching fraction along with its CP asymmetry. Since 
the branching fractions and CP asymmetries of all the 
other pn final states have been measured @], p°7r° is the 
only channel that remains to complete the isospin pen- 



tagons. A simplification, whereby the pentagons collapse 
into quadrangles, is also possible if the p°Tr° amplitude is 
sufficiently small. 

An alternative technique to measure fa from the pit 
system, even if penguin contamination is large, is a time- 
dependent amplitude analysis of B° — > 7r + 7r~7r° 0. 
Here, the interferences between p + ir~ , p°ir° and p~ir + 
contributions to the 7r + 7r~7r° final state provide the crit- 
ical information on the unknown phases introduced by 
penguin amplitudes. Recently, the first time-dependent 
studies of the ir + ir~ir° Dalitz plot have been per- 
formed 0]. In these studies, a simplification is made with 
the assumption that the p°ir° contribution is small. A 
more complex time-dependent Dalitz analysis is required 
if this is not the case. 

The Belle Collaboration reported first evidence of the 
B° — > p°7r° decay Q with a branching fraction larger 
than most predictions 0, and a central value above the 
90% confidence-level upper limit set by the BaBar Col- 
laboration jlOj. In this paper, we report an improved 
measurement of the B° — > p°ir° branching fraction 
using 2.5 times more data, and perform a first direct CP 
violation search in this mode. The results are consistent 
with and supersede those reported in our previous pub- 
lication. The analysis is based on (385.8 ± 4.8) x 10 6 
BB pairs, collected with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider ^3] that operates at the 
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T(4S f ) resonance. The production rates of B + B and 
B°B° pairs are assumed to be equal. 

The Belle detector 0, Q] is a large-solid-angle mag- 
netic spectrometer that consists of a silicon vertex detec- 
tor, a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of 
CsI(Tl) crystals (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to 
detect K Q L mesons and to identify muons. 

The B reconstruction procedure is identical to our pre- 
viously published analysis @- Charged tracks are re- 
quired to originate from the interaction point and have 
transverse momenta greater than 100 MeV/c. Pions 
are identified by combining information from the ACC, 
TOF and the CDC dE/dx measurements. We further 
reject tracks that are consistent with an electron hypoth- 
esis. Pairs of photons with invariant masses in the range 
0.115 GeV/c 2 < to 77 < 0.154 GeV/c 2 are used to form 
the 7T° mesons. The photon energy in the laboratory 
frame is required to be greater than 50 (100) MeV in 
the barrel (endcap) region of the ECL. The tt° candi- 
dates are required to have transverse momenta greater 
than 100 MeV/c in the laboratory frame and a loose re- 
quirement is made on x 2 , the goodness of fit of a n° 
mass-constrained fit of the two photons. We also veto 
possible contributions to 7r + 7r~7r from charmed (b — > c) 
decays: B° -> D~n + , D°tt and J/ipir°. 

Signal B candidates are identified with two kinematic 
variables: the beam-energy constrained mass M bc = 
V ^bcam/ c4 P% l c2 an d tne ener gy difference AE = 
Eb — -Ebcam- Here, Eb (pb) is the reconstructed energy 
(momentum) of the B candidate, and E^ eam is the beam 
energy, all expressed in the centre-of-mass (CM) frame. 
We consider candidate events in the region —0.2 GeV < 
AE < 0.4 GeV and M bc > 5.23 GeV/c 2 ; and define sig- 
nal regions in AE and M bc as -0.135 GeV < AE < 
0.082 GeV and 5.269 GeV/c 2 < M bc < 5.290 GeV/c 2 . 
To select p°ir° from the ir + Tr~ir° candidates, we require 
the 7r + 7r~ invariant mass to be in the range 0.5 GeV/c 2 < 
mjr+Tr- < 1.1 GeV/c 2 and the p° helicity angle to sat- 
isfy | 

cos ^hcil ^ 0-5; where 0^ el is defined as the angle 
between the negative pion direction and the opposite 
of the B direction in the p rest frame. We explicitly 
veto contributions from B a — > p^n^ by the require- 
ment m w ± 7T o > 1.1 GeV/c 2 . This requirement also ve- 
toes the region of the Dalitz plot where the interference 
between p°tt° and p^n^ is strongest. After all selection 
requirements, 11% of events have more than one candi- 
date. Among those candidates the one with the smallest 
xltx/ nd f + xlo I ndf is selected, where Xvtx is tne good- 
ness of fit of a vertex-constrained fit of 7r + 7r~ . 

The dominant background originates from continuum 
e + e~ —> qq (q = u, d, s, c) production. To separate the 



jet-like qq events, we use event shape variables: five mod- 
ified Fox- Wolfram moments 0] , combined into a Fisher 
discriminant. We further combine the cosine of the B 
meson flight direction in the CM system with the out- 
put of the Fisher discriminant into a signal/background 
likelihood variable, £ s /&, and define the likelihood ratio 
1Z = C s / (C s + Cb). Additional discrimination against 
continuum is achieved through use of the 6-flavour tag- 
ging algorithm We use the parameter r, with values 
between and 1, as a measure of the confidence that the 
remaining particles in the event (other than 7r + 7r~7r°) 
originated from a flavour specific B meson decay and — 
as a corollary — not from a continuum process. 

We use an iterative procedure to find the optimal con- 
tiguous area in r-lZ space by maximising N s /^/N s + Nb, 
where N s (Nb) is the expected number of signal (back- 
ground) events in the AE and M bc signal regions. Here, 
the optimisation procedure assumes a branching fraction 
for B° -> p 7r of 3.3 x 10~ 6 Anticipating the use 

of r for its primary purpose of flavour tagging in CP 
asymmetry fits, the borders of the contiguous area were 
constrained to match the six r bins employed in previous 
analyses. The result of the optimisation procedure is that 
we select events within the region shown in Fig.^a). We 
find 1397 candidates remain in the data. 

We obtain the signal yield using an extended unbinned 
maximum-likelihood fit to the AE-M^ C distribution of 
the selected candidate events. The likelihood function is 
defined as 

£ = ex P (- £ N jt )j J] \T Nj,^ . (1) 

Here, the index i is the event identifier; I distinguishes 
events in various r bins; and j runs over all six compo- 
nents included in the fitting function — one for the signal, 
and the others for continuum, b — > c combinatorial, and 
the charmless B backgrounds: B + — > p + p a , B + — > p + ir° 
and B + — * 7r + 7r°. Njj, represents the number of events, 
and Vj — Pj(M^ c , AE 1 ) are two-dimensional probability 
density functions (PDFs). 

The PDFs for signal, b — + c and charmless B back- 
grounds are taken from smoothed two-dimensional his- 
tograms obtained from Monte Carlo (MC) simulations. 
For the B + — > p + p° channel, we assume a 100% longi- 
tudinally polarised decay ^^|. Small corrections to MC 
peak positions and widths are applied to the signal PDF. 
These factors are derived from control samples of recon- 
structed decays B° -> D*~ p+ (D*~ -> ~D a ir- , D° -> 
K+tt-; p+ -> tt+tt ) and B+ -> D°p+ (D° -> A'+tt"; 
p + — * 7T + 7r°), in which we require that the ir° momentum 
be greater than 1.8 GeV/c in order to mimic the high mo- 
mentum 7T° in our signal. The two-dimensional PDF for 
the continuum background is described as the product 
of a first-order polynomial in AE and an ARGUS func- 
tion [l9j in M bc . All of the shape parameters describing 
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FIG. 1: (a) Distribution of signal (continuum) events in r-lZ space shown with open (shaded) proportional boxes; the marked 
region (top-right) indicates the area selected, (b) (c) Distribution of AE(Mbc) in the signal region of A/bc(A£'). Projection 
of the fit result is shown as the thick solid curve; the thin solid line represents the signal component; the dashed, dotted and 
dash-dotted curves represent, respectively, the cumulative background components from continuum processes, b — » c decays, 
and charmless B backgrounds, (d), (e) Distributions of fit yields in m n + 7T - and cos#£ cl variables for p°7r° candidate events. 
Points with error bars represent data fit results, and the histograms show signal MC expectation; the selection requirements 
described in the text are shown as dashed lines. 



the continuum background are free parameters in the fit. 
The normalisations of B + — > p + Tr° (21.7±4.4 events) and 
B + — > 7T + 7r° (2I.0±5.5 events) are fixed in the fit accord- 
ing to previous measurements 0, l2(l | , and that of b — > c 
background (62 ± 62 events) according to MC expecta- 
tion (assigning a conservative error); the normalisations 
of all other components arc allowed to float. 

The fit result is shown in Fig. ^b) and (c). The 
signal yield is found to be 50.91^3 4 with 4.5er signifi- 
cance. The significance is defined as ^J— 2 ln(£o/£max), 
where £ max (£q) denotes the likelihood with the signal 
yield at its nominal value (fixed to zero). The contri- 
bution from B + — > p + p° decays (which peaks in the 
low AE region) is obtained from the fit as 43.1^2' 2 
events; this value is consistent with the MC expectation 
•9ig'g events) based on our branching fraction mea- 
surement of B + — > p + p° 01- A possible contribution 
from B — > w(7r + 7r~7r°)7r° decays is also accounted for 
by floating the B + — > p + p° PDF, since the two decays 
have similar distributions in AE and M^ c . To verify that 
the signal candidates originate from B° — > p°ir° decays, 
we change the criteria on m^+Tr- and cos6^ cl in turn, 
and repeat fits to the AE-M-^ C distribution. The yields 
obtained in each 771^+^- and cos(9^ el bin are shown in 
Fig-dd) and (e). 

The cos 6^ e [ distribution is used to limit contributions 
from 5° -> an , / (980)7r°, rj' n° , K s tt° and vr+ t^tt 
(nonresonant), which are expected to be flat in this vari- 
able. We perform a \ 2 fit including components for pseu- 
doscalar — * pseudoscalar vector (PV ~ cos 2 #£ cl ), and 
pseudoscalar — > pseudoscalar scalar (PS ~ flat) decays, 
for which the shapes are obtained from our p°n° signal 
MC, and a sample of an MC |2l|, respectively. We also 
include a linear term to allow for possible interference. 
We find that the PS level is consistent with zero; tak- 
ing its uncertainty into account, we assign a model error 
of ii5°o% to the PV component. The m n + n - distribu- 
tion is consistent with the expectation from B — > p°n° 



production. 

To extract the branching fraction, we determine the 
reconstruction efficiency, (4.99 ± 0.03)%, from MC and 
correct for small differences between data and MC in 
the pion identification and continuum suppression re- 
quirements. The correction factor due to charged pion 
identification (0.872) is obtained in bins of track momen- 
tum and polar angle from an inclusive D* control sam- 
ple (D*~ — ► D°n~, D° — > K + n~). The corresponding 
systematic error is ±3.1%. For the continuum suppres- 
sion requirement on r and 1Z, we use the control sample 
B° -> D~p+ (£)- -> K+ir-ir-; p+ -> tt+tt ) to obtain 
an efficiency correction factor of 0.972 and a correspond- 
ing systematic error of ±6.0%. 

We calculate additional systematic errors from the fol- 
lowing sources: PDF shapes by varying parameters by 
±Icr (+2 '0%); 7T° reconstruction efficiency by comparing 
the yields of r\ — > tt°tt°tt and r\ — ► 77 between data 
and MC (±4.0%); track finding efficiency from a study 
of partially reconstructed D* decays (±2.4%); and data- 
MC efficiency differences due to the AE > — 0.2 GeV re- 
quirement (±2.0%). We repeat the fit after changing the 
normalisation of the fixed background components ac- 
cording to the given errors and obtain a systematic error 
of iYg%- Using a large MC sample, the total systematic 
error from possible charmless B decays not otherwise in- 
cluded, B° -► K*°n° (5.4%), B+ -> K*+n° (1.5%) and 
B° -> K+p- (0.5%), is ±5.6%. 

When the normalisations of all the backgrounds fixed 
in the fit are simultaneously increased by la, the sta- 
tistical significance decreases from 4.5<r to 4.2er; we in- 
terpret the latter value as the significance of our re- 
sult. Finally, we estimate the uncertainty due to possi- 
ble interference with B° — > p^n^ by varying the m,± I o 
veto requirement from m 7T ± 7r o > MeV/c 2 (no veto) to 
m^i^o > 1.7 GeV/c 2 . We find the largest change in the 
result to be within ±16%, and we include this value in the 
model error, so that the obtained B° — > p°n° branching 
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fraction is 

B= (3.12±g;||(stat) ± 0.33(syst)±g;j^ (model)) x 10~ 6 . 

Having observed a significant B — > p°ir a signal, we 
utilize the B° / B° separation provided by the flavour tag- 
ging to measure the CP asymmetry. For this purpose we 
replace Vj of Eq. (1) with the expression 



1 + q l ■ (A 



CP)j,l 



P^M^AE* 



(2) 



in which the indices keep the same meaning. In this equa- 
tion, q represents the b- flavour charge [q = +1 (—1) when 
the tagging B meson is a B° (B )] and A CP denotes 
the effective charge asymmetry, such that {A CP )j y i = 
(Acp)j(l - 2xd)(l - 2wi). Here, (A C p)j are the charge 
asymmetries for the signal and the ba ckgr ound compo- 
nents. Further, x<i = 0.182 ± 0.015 Q is the time- 
integrated mixing parameter and wi is the wrong-tag 
fraction. For continuum background, \d and wi are set 
to zero. The data is divided into the six r-bins, and the 
r-dependent wrong-tag fractions, wi (I = 1,...,6), are 
determined using a high statistics sample of self-tagged 
B° -> £>M-7r+, D*-p+ and D*-£+v events 0. 

The total number of signal, continuum background and 
p + p a events are free parameters in the fit, and the re- 
maining background components (from b — * c, p + ir° and 
7r + 7r decays) are fixed. Also, the relative fractions for 
the signal and continuum background components in dif- 
ferent r bins are allowed to float in the fit; for the b — > c 
and charmless B decay backgrounds, they are fixed. The 
only free Acp parameter in the nominal fit is that of 
our signal; the others are fixed to be zero (for continuum 
and b — > c) or at their previously measured values (for 
charmless B backgrounds) j^Jj. We measure the direct 
CP asymmetry in B° — > p°ir° decays to be 



A 



CP 



0.53t^(stat)^^(syst). 



The impact of background asymmetry (1° 127) is the 
largest contribution to the systematic error; it is esti- 
mated by releasing, in turn, all of the background Acp 
parameters (limiting them within ±1ct range of their 
measured values for the charmless B decays), and sum- 
ming in quadrature the differences obtained from the cen- 
tral Acp value. A similar sum gives igioi? as * ne s y s ~ 
tematic uncertainty obtained by varying all other fixed 
parameters in the fit, including \d and u>i values, by ±ler. 
Finally a systematic error of ±0.058 is obtained as a re- 
sult of a null asymmetry test, when the same analysis pro- 
cedure is applied to the B° — > D~ p + (D~ — » K + n~ir~; 
p + — > 7T + 7r°) control sample. To illustrate the asymme- 
try, we show the results separately for p°7r° candidate 
events tagged as q = +1 and q = — 1 in Fig.|3 

In summary, using 386 x 10 6 BB pairs, we confirm 
evidence of B a — > p°ir° decays with a branching fraction 
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FIG. 2: AE and M^ c distributions (with projections of the fit 
results) shown separately for events tagged as q = +1 (left) 
and q = — 1 (right). 



higher than most theoretical predictions [9| . The central 
value remains only slightly above the 90% confidence- 
level upper limit set by the BaBar Collaboration [Io| . 
and is in agreement with the upper limit set by the CLEO 
Collaboration 0. Our measurement is consistent with, 
and supersedes, our previous result [8j. We have also 
performed a first measurement of direct CP violation in 
the B° — > p°7T mode and find no statistically significant 
asymmetry. 

The large p°ir° branching fraction suggest that one 
can only impose a loose constraint on penguin uncer- 
tainty in the determination of <f>2 from time-dependent 
B° — > p ± tt t measurements. It also implies that a useful 
measurement of <f>2 from the full pir isospin analysis may 
be impractical even with super B-factory like luminosi- 
ties |23j. Therefore, we can expect that the best measure- 
ments of <f>2 from the pir system will come from the full 



time-dependent amplitude analysis of B° — > 7r + 7r~7r°. 

We thank the KEKB group for excellent operation of 
the accelerator, the KEK cryogenics group for efficient 
solenoid operations, and the KEK computer group and 
the Nil for valuable computing and Super-SINET net- 
work support. We acknowledge support from MEXT and 
JSPS (Japan); ARC and DEST (Australia); NSFC and 
KIP of CAS (contract No. 10575109 and IHEP-U-503, 
China); DST (India); the BK21 program of MOEHRD, 
and the CHEP SRC and BR (grant No. R01-2005-000- 
10089-0) programs of KOSEF (Korea); KBN (contract 
No. 2P03B 01324, Poland); MIST (Russia); MHEST 
(Slovenia); SNSF (Switzerland); NSC and MOE (Tai- 



6 



wan); and DOE (USA). 



[1] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 
652 (1973). 

[2] K. Abe et al. (Belle Collab.), Phys. Rev. Lett. 95, 101801 
(2005); B. Aubert et al. (BaBar Collab.), Phys. Rev. 
Lett. 95, 151803 (2 005); A. Somov et al. (Belle Col- 
lab.), hep-ex/0601024 to appear in Phys. Rev. Lett; 
B. Aubert et al. (BaBar Collab.). Phys. Rev. Lett. 95, 
041805 (2005). 

[3] C. C. Wang et al. (Belle Collab.), Phys. Rev. Lett. 94, 
121801 (2005). 

[4] M. Gronau and D. London, Phys. Rev. Lett. 65, 3381 
(1990); H. J. Lipkin, Y. Nir, H. R. Quinn and A. E. Sny- 
der, Phys. Rev. D 44, 1454 (1991); Y. Grossman and H. 
R. Quinn, Phys. Rev. D 58, 017504 (1998). M. Gronau, 
D. London, N. Sinha and R. Sinha, Phys. Lett. B 514, 
315 (2001). 

[5] A. E. Snyder and H. R. Quinn, Phys. Rev. D 48, 2139 
(1993). 

[6] A. Gordon et al. (Belle Collab.), Phys. Lett. B 542, 183 
(2002); B. Aubert et al. (BaBar Collab.), Phys. Rev. Lett. 
91, 201802 (2003); J. Zhang et al. (Belle Collab.), Phys. 
Rev. Lett. 94, 031801 (2005); C.-P. Jessop et al. (CLEO 
Collab.), Phys. Rev. Lett. 85, 2881 (2000). 

[7] B. Aubert et al. (BaBar Collab.), hep-ex/0408099 sub- 
mitted to the ICHEP 04 conference. 

[8] J. Dragic et al. (Belle Collab.), Phys. Rev. Lett. 93, 
131802 (2004). 



[9] C.-W. Chiang, M. Gronau, Z. Luo, J. L. Rosner and D. 
A. Suprun, Phys. Rev. D 69, 034001 (2004); M. Beneke 
and M. Neubert, Nucl. Phys. B 675, 333 (2003). 
[10] B. Aubert et al. (BaBar Collab.), Phys. Rev. Lett. 93, 
051802 (2004). 

[11] Throughout the paper, the inclusion of the charge con- 
jugate mode decay is implied unless otherwise stated. 

[12] S. Kurokawa and E. Kikutani, Nucl. Instr. and Meth. A 
499, 1 (2003), and other papers included in this volume. 

[13] A. Abashian et al. (Belle Collab.), Nucl. Instr. and Meth. 
A 479, 117 (2002). 

[14] Y. Ushiroda (Belle SVD2 Group), Nucl. Instr. and Meth. 
A 511, 6 (2003). 

[15] G. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 

[16] H. Kakuno et al, Nucl. Instr. and Meth. A 533, 516 
(2004). 

[17] The number is the average of central values reported in 

references [il Ho|. 
[18] J. Zhang et al. (Belle Collab.), Phys. Rev. Lett. 91, 

221801 (2003). 

[19] H. Albrecht et al. (ARGUS Collab.), Phys. Lett. B 241, 
278 (1990). 

[20] S. Eidelman et al. (Particle Data Group), Phys. Lett. B 

592, 1 (2004) and 2005 partial update for edition 2006. 
[21] We model the a particle with a mass of 478 MeV/c 2 and 

a width of 324 MeV/c 2 , based on E. M. Aitala et al. 

(E791 Collab.), Phys. Rev. Lett. 86, 770 (2001). 
[22] Heavy Flavor Averaging Group (HFAG), "Averages of b- 

hadron Properties as of Winter 2005", hep-ex/0505100 
[23] J. Stark, eC onf C0304052, WG423 (2003), 

hep-ph/0307032 



